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ABSTRACT 


The  basic  objective  of  the  work  reported  herein  is  to  study,  describe 
and  measure  the  transport  of  nuclear  radiation  through  air  ducts  having  walls 
of  different  materials.  Particular  emphasis  is  placed  on  full  scale  concrete 
ducts  yrtiich  are  intended  for  use  as  underground  shelter  entranceways  and  contain 
a  right  angle  bend. 

e 

It  is  shown  that  a  modified  albedo  theory  will  rather  accurately 
describe  the  transmission  of  both  gammas  and  neutrons  in  lead  and  concrete 

t 

ducts.  Theory  is  con^jared  with  experimental  measurements  of  gamma  dose  and 
neutron  flux  attenuation  factors.  Neutron  number  albedo  measurements  are  made 
and  reported. 

Work  by  various  theorists  and  experimenters  are  included  for  complete¬ 
ness  and  comparison.  This  program  rather  clearly  points  out  direction  future 
programs  should  take  and  the  areas  in  which  major  emphasis  is  needed.  ^ 
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MDUnON  STBEHnm  IN  SHELTER  ENIBANCENAYS 


I.  INTRODUCTION 

In  the  winter  of  1959,  the  Armour  Research  Foundation  began  work  on  a 
program  designed  to  develop  basic  Information  on  the  transport  of  nuclear 
radiation  through  ducts  and  shelter  entranceways.  The  need  for  such  a 
program  is  clearly  In  evidence  Arom  the  shelter  radiation  level  measure¬ 
ments  made  during  the  Teapot  and  Pluabob  test  series.  The  work  reported 
herein  Is  a  combined  theoretical  and  experimental  program  designed  to 
Investigate  both  neutron  and  gama  ray  transport  characteristics. 

HVhlle  the  basic  objective  was  information  necessary  for  adequate 
radiation  protection  design  of  personnel  shelters,  it  was  recognised  that 
because  no  experimental  data  existed  for  full  scale  entranceways  a  eontrlbu>- 
tion  to  the  field  of  radiation  transport  could  be  made.  Further,  the  absence 
of  experimental  work,  In  general,  adds  value  to  almost  any  measurements 
which  are  made.  Various  authors  have  reported  analytical  receipts  vhleh 
invoke  the  albedo  concept.  Numerous  calculations  of  albedos  have  been 
reported  but  the  paucity  of  experimental  values  renders  comparison  of  theoiy 
and  experiment  virtually  Impossible  In  all  but  a  few  special  casts. 

Consequently,  this  program,  though  not  part  of  the  original  planning, 
attempts  to  determine  experimental  albedo  values.  We  were  successful  only 
with  neutron  number  albedo,  other  atten?5t5  being  frustrated  by  scattering  from 
undesired  surfaces,  inadequate  sources  or  detectors,  etc. 


This  report,  in  addition  to  describing  the  work  done  on  this  program, 
attempts  to  pull  together  most  work  done  to  date  of  which  we  have  knowledge. 


II.  SIIIMARY  OF  PREVIOUS  ffQBK 
A.  Analytical 


1,  Simon  and  Clifford^ 

Simon  and  Clifford  presont  the  following  expression  for  the 
uncollided  fl\ix  at  the  end  of  a  straight  cylindrical  duet  from  a  plane  source 
of  strength 


lAiere 


S 


radius 


length  of  duet 

2 

strength  in  n/cm  -  sac  isotropically  In  the 
forward  direction. 


This  formula  ia  valid  for  ducts  where  the  radius  S  is  small  compared  to 
the  length  ,  The  scattered  flux  from  reflections  from  the  walls  is  given 


by: 


N. 


s 


2mF 


(2) 


^Slmonj  A.  and  Clifford,  C.  E.  ’’The  Attenuation  of  Neutrons  by  Air  Ducts  in 
Shields.”  Nuc.  Sci.  and  Engr.  1,  156-166  (1956). 


YAiere 


agd  X  is  the  reflection  coefficient  (albedo). 


The  dose  albedo  is  defined  as  the  ratio  of  the  dose  from  a 
surface  to  the  dose  impinging  upon  the  surface.  If  9^  is  independent  of 
angle,  equation  (2)  reduces  to  j 


Simon  and  Clifford  at  this  point  enforces  their  assumption  of 
approximates  the  integral  as  from  0  to  l/2  instead  of  from  0  to  1. 
technique,  the  scattered  flux  reduces  to 


1  and 
By  this 


(h) 


If  the  albedo  is  part  isotropic  (A)  and  part  cosine  (B)  the  total  flux  becomes 


1  B 


(5) 


It  might  be  noted  that  the  Integral  in  equation  (3)  can  be 
2 

integrated  directly J  however,  the  solution  is  very  long  and  since  the 
assumption  that  the  albedo  does  not  change  with  angle  is  poor,  the  added 
value  of  the  direct  Integration  is  highly  doubtfhl.  The  solution  shown  in 
equation  (1;)  is  probably  as  good  as  is  Justified  by  the  other  assunptions 
used  in  the  derivation. 

The  dose  after  a  bend  is  equivalent  to  two  ducts  Joined 


at  an  angle  of  0.  The  equation  for  the  dose  is 


where  k  is  a  factor  converting  flux  to  dose.  The  use  of  an  angular  dependent 


albedo  in  equation  (6)  increases  its  accuracy  but  the  assumption  of^^^l 
limits  its  usefulness. 

Simon  and  Clifford  also  present  a  single  scattering  approach 
for  the  attenuation  of  neutrons  down  a  straight  duct.  They  assume  isotropic 
scattering  and  that  the  number  of  neutrons  at  a  point  located  at  a  distance  x 
into  the  wall  is  the  number  incident  on  the  surface  multiplied  by  exp  (- 
Their  expressions  for  the  scattered  flux  at  a  point  inside  a  straight 
cylindrical  duct  is  :  -f  ^  ^  j 

Cb_  tils  I  Mx->-S)e  * 

’  W  /  fttU  i*xt  1  i  +xf'  (7) 


2 

ARF  Monthly  Report  Number  6  of  this  contract- Integration  performed  by  H.  Wilf. 


irtiere 


and  X  ■  depth  into  surface,  y  •  distance  from  source, is  scattering  cross 
section  and  is  the  total  cross  section.  Ihe  geometry  of  the  problem  is 

X 

shown  in  lig.  1. 


K - / 


Fig.  1.  Single  Scattering  Geometry 


Since  the  function  under  the  Integral  in  equation  (7)  peaks  at  x  -  0,  all 
slowly  varyflng  functions  are  evaluated  at  x  <■  0.  Setting  yK<  1  and 


1  we  obtain  for  the  scattered  flux 


Ts 

t: 


(9) 


Ihis  shows  that  the  scattered  flux  varies  as  the  cube  of  the  length,  not  as 
the  square  as  was  indicated  from  the  albedo  method  (see  equation  (ii)).  It  can 
be  shown  that  equations  (It)  and  (9)  converge  when  the  angular  dependant  albedo 
is  used  in  equation  (U).  The  agreement  is  good  for  long  leg  lengths. 

Slsion  and  Clifford  present  also  an  approximete  derivation  for 
multiple  reflection  which  can  be  summarised  in  one  statement  thusly:  for  a 
point  surrounded  by  reflecting  surfaces,  replace  the  albtdo  wlth-^^.  This 
statement  involves  a  good  many  assumptions  and  approKimations  but  does 
yield  good  results  for  snisll  ^  and  fair  iccurecy  for  ^  up  to/v'l/S.  The 
single  scattering  technique  at  outlined  by  Simon  and  Clifford  will  not  work 
very  well  for  gamma  rays  because  of  tha  large  energy  degradation  and  high 
forward  scattering. 


-6- 


2,  Roe^ 

Roe  treats  the  streaming  problem  liy  a  one-group  diffusion 
approach.  He  uses  the  same  limiting  assumption  that  the  radius  is  small 
compared  to  the  leg  length.  Roe  does  reach  a  fairly  rigorous  series  solution 
for  very  large  leg  lengths.  His  results  for  the  scattered  flux  are* 

-  (XO) 

where 

/^tK 

(K.S) 


/ 

=  leg  length 
5  =  duct  radius 

K  and  IL,  are  the  standard  Bessel  functions. 

0  1  ^ 

The  results  arv  fairly  good  for  large  /y  but  suffers  from  the  usual  diffusion 
theory  limitations  near  boundaries. 


3 

Roe,  G.  M.  "The  Penetration  of  Neutrons  Through  an  Empty  Qylindrical  Duct 
in  a  Shield."  KAPIr-712  (1952). 
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3.  Barcua^ 

Barcus  starts  In  the  same  manner  as  Simon  and  Clifford, 
hovrever  he  uses  an  Integral  approach  to  the  problem.  He  uses  the  assump¬ 
tion  of  Simon  and  Clifford  in  that  the  radius  is  very  small  as  coti^ared  to 
the  leg  length.  The  main  difference  is  that  the  albedo  concept  Is  not  used. 

Ihe  expression  for  the  direct  flux  at  point  x,  Triiich  is  a 
distance  Jf  away  from  a  plain  isotropic  sourc^is 


This  expression  reduces  to  the  Simon  and  Clifford  results  if  we  let  be 
much  less  than  one.  Barcus  then  derives  an  expression  for  the  scattered 
component.  He  used  a  Kernel  approach  with  the  transport  equation.  His 
final  expression  for  the  scattered  component  is: 


^  Barcus,  J.  R.  "Transmission  of  Neutrons  by  Cylindrical  Ducts  Penetrating 
Radiation  Shields."  TID-ii$00  (lUth  Ed.)  Physics  and  Mathematics. 
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vihere 


}CXstd^^i  j  ^  ,  O^yLoi  ^  *  mean  free  path 


For  a  90®  band,  equation  (11)  reduces  to 


5?= 


(13) 

for  the  scattered  flux. 
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d 

k,  LeDoux  and  Chilton 


This  work  was  received  as  this  report  was  in  preparation 
and  as  such,  no  time  was  available  for  complete  analysis  or  checking  with 
the  work  of  others.  The  albedo  approach  is  used  with  the  usual  limitation 


of 


Iheir  results  appear  in  the  form  of  four  contributing  terms* 

(a)  “  radiation  scattered  from  areas  A^,  A^, 
and  A|^  -  utilizing  albedo  reflection.  These  areas 
can  see  both  source  and  detector. 

(b)  -  radiation  scattered  from  areas  A^,  A^,  A^,  and 
‘  Ag  >  utilizing  albedo  reflection.  These  are  areas 

idiere  the  radiation  has  passed  through  the  corner, 

(c)  0  -  corner  scattering  of  radiation  \\diich  is  then 

s 

.  directed  toward  detector, 

(d)  D.  -  direct  corner  penetration  from  source  to  detector. 
The  rule  for  inclusion  in  the  solution  of  the  various  terms 

is  based  on  a  greater  than  10  °/o  effect.  Such  a  cutoff  point  is  taken  in 
order  to  limit  the  complexity  of  the  solution, 

,  The  summations  are  performed  using  average  distances  and 
average  angles  to  determine  the  albedo.  Tables  are  utilized  to  determine  the 
attenuation  ratios  for  various  rectangular  geometries. 


^  LeDoux,  J.  C.  and  Chilton,  A.  B,  "Attenuation  of  Gamma  Radiation  Through 
Two-legged  Rectangular  Dirts  and  Shelter  Entranceways  -  An  Analytical 
Approach.”  Soon  to  be  released  under  the  auspices  of  the  U.  S.  Naval  Civil 
Engineering  Laboratory,  Port  Heuneme,  California, 


-10- 


B.  Experimental 

As  a  convenient  reference  to  the  reader,  we  have  drawn 
together  that  experimental  data  of  udilch  we  are  aware.  The  most  complete  work 
done  Is  by  Elsenhauer^  at  the  National  Bureau  of  Standards  In  which  gamma 
attenuation  measiu*ements  are  made  In  small  concrete  ducts.  He  Investigated 
the  "corner  effect"  by  substituting  lead  for  concrete.  Unfortunately,  his 
measurements  are  all  made  with  pocket  dosimeters  with  attendant  low  accuracy 
and  reproducibility.  Figure  2  shows  the  geometry  and  his  results  are  presented 
in  Table  I.  We  note  that  the  corner  effect  is  more  pronounced  in  the  smaller 
duct  and  effects  points  near  the  comer  most.  When  a  leg  length  (D)  of  68. 7  cm 
is  reached  the  corner  effect  still  contributes  32  °/o  (1.97  mr/hr  vs  1,33). 

This  indicates  that  for  small  ducts,  leakage  through  and  scattering  in,  the 
corner  has  a  significant  effect  on  what  arrives  at  the  exit  of  duct.  It  is 

shown  that  this  same  conclusion  holds  for  full  scale  entranceways. 

7 

Hungerford  has  collected  a  series  of  experimental  albedo 
measurements  which  attempt  to  obtain  a  total  albedo  for  gamma  rays  and  neutrons. 
Reflection  from  a  flat  surface  is  measured  by  a  detector  shielded  from  the 
direct  source  radiation,  Hungerford  reports  values  of  0.12  for  PoBe  neutrons 
and  0. OU  for  Co-60  gamma  rays  in  concrete. 


Private  communication:  Eisenhauer  to  LeDoux. 

7 

Hungerford,  H.  E.  "Some  Scattering  Experiments  Performed  at  the  Bulk 
Shieldinr  Facility."  C?-$2^k-99 ,  April.  16,  1952  (OTS). 
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8  9 

Work  by  Cure  and  Hurst  and  by  Strickler,  Gilbert  and  Auxier 

using  techniques  similar  to  Hungerford  report  a  dose  albedo  value  of  0.32  for 

PoBe  neutrons. 


III. 


ANALYTICAL  PROGRAM 


As  described  earlier,  a  number  of  theoretical  recipes  have  been 
developed  to  describe  the  transport  of  nuclear  radiation  through  straight 
diicts  and  through  ducts  having  bends.  All  of  these  previous  recipes  suffer 
from  inherent  approximation  which  seriously  limits  their  applicability.  The 
approximations  are  in  the  form  of  geometry  limitations  (i.  e. , 
energy  independence,  angular  independence,  use  of  diffusion  theory,  single 
scattering  approximation  (dropping  of  higher  order  scattering)  and  others. 

Experience  with  similar  problems  in  reactor  theory  predicts  that 
use  of  integral  transport  theory  in  which  angle  and  energy  dependent  scattering 
is  employed  would  provide  a  superior  recipe.  Unfortunately,  so  many  unknowns 
exist  that  such  an  elaborate  approach  cannot  be  justified  at  this  time. 
Similarly,  a  Monte  Carlo  approach  is  not  justifiable  under  this  program. 

Our  approach  has  been  to  investigate  the  more  simple  treatments 
with  the  goal  of  enhancing  their  value  by  suitable  correction  terms,  much  in 
the  manner  that  diffusion  theory  may  be  improved  by  making  certain  transport 
corrections. 


Cure,  J,  W.  and  Hurst,  G,  S.  Nucleonics  12.  36-38 

Strickler,  T.  D. ,  Gilbert,  H.  E.  and  Auxier,  J,  A.  "Fast  Neutron  Scattering 
from  Slabs,  "  Nuc.  Sci.  Engrg.  3.  11-18  (1957). 


-Ill- 


A.  Aiialytic_al  Methods 

The  first  method  investigated  Is  single  scattering.  The  assumption 

Is  made  that  the  dose  at  a  point  is  due  to  a  direct  imcoUided  flux  (D^)  and 

to  a  flux  once  scattered  (D  )  from  the  walls.  The  assumption  of  single 

8 

scattering  In  the  walls  Implies  that  only  the  once  scattered  flux  makes  a 
significant  contribution  to  the  dose.  In  general  vrtxen  ^ 

sufficient  condition  for  the  single  scattering  approximation. 

If  this  is  not  satisfied  one  might  expect  that  the  amount  of 
radiation  predicted  by  single  scattering  will  be  an  underestimate  and  this  was 
confirmed  by  experiment.  We  therefore  have  to  include  a  correction  due  to 
multiple  scattering.  Even  double  scattering  would  be  very  laborious,  for  all 
but  the  simplest  cases,  hence,  we  look  for  a  more  phenomenological  correction 
to  single  scattering.  One  way  of  doing  so  is  to  invoke  a  buildup  factor.  This 
is  a  number  relating  the  amount  of  radiation  arriving  at  a  certain  point,  .to 
that  expected  if  no  scattering  occurred.  The  buildup  factor  is  a  function  of 
the  number  of  mean  free  paths  the  radiation  has  traveled  in  the  scattering 
material  and  has  the  effect  of  increasing  the  amount  of  radiation  reaching 
the  scattering  point. 

An  equivalent  method  in  the  case  of  neutrons  is  obtained  by  the 
introduction  of  the  removal  cross  section.  This  is  defined  as  a  measure  of 
the  ability  of  the  material  to  remove  neutrons  from  a  beam. 

Another  approach  is  known  as  the  albedo  concept.  The  dose  albedo 
is  defined  as  the  ratio  of  the  dose  emerging  from  a  certain  point  on  the 
surface,  to  the  dose  falling  on  the  surface.  The  dose  and  number  albedos 
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are  calculated  for  gamma  rays  as  a  function  of  energy  and  angle  of  incidence, 

> 

Unfortunately  the  number  albedo  for  neutrons,  which  is  necessary  for  the 
type  of  neutron  detector  we  emploj^  is  not  available.  For  this  reason  several 
experiments  are  performed  and  an  analysis  of  the  data  yields  a  consistent 
albedo,  vrtiich  is  assumed  to  be  constant  over  the  surface  as  well  as  over  energy. 
1.  Single  Scattering 

We  will  now  present  a  calculation  based  on  single  scattering. 
We  wish  to  find  the  radiation  at  a  detector  due  to  a  source  S;  and  we  will 
consider  two  typical  detector  pointss  a)  points  which  can  see  the  source  (D'), 
and  b)  points  which  get  only  scattered  radiation  (D*’). 


Fig.  3.  Single  Scattering  Geometry 


Figure  3  shows  a  top  view  of  the  entrancewsiy-.  The  contribu¬ 


tion  from  a  point  source  to  point  D'  will  be 

D  — 

Here  5  measures  the  strength  of  the  isotropic  point  source  and  k  is  a  con¬ 
version  factor  from  flux  to  dose. 

The  contribution  to  the  dose  at  D',  from  an  element  of  volume 
dV  in  the  wall,  will  be; 


Ci5) 


Here  2^  is  the  total  cross  section  for  the  source  energy,  ^2^  the  same 

quantity  for  a  flux  with  degraded  energy.  ^  is  the  macroscopic  scattering 

s 

cross  section  which  is  not  necessarily  angular  independent.  The  part  in  the 
brackets  of  (l5)  represents  the  source  per  unit  volume  of  the  once  scattered 
radiation.  To  obtain  the  total  contribution  due  to  single  scattering  one  has 
to  integrate  over  that  part  of  the  wall  volume  which  can  be  seen  by  the  de¬ 
tector.  The  scattered  contribution  (D  )(2'  is  angularly  independent)  becomes 

s  s 


-17- 


The  volume  integration  may  be  reduced  to  a  surface  integration 
by  approximation.  From  simple  geometrical  considerations  one  obta'ins 


(17) 


(18) 


The  exponential  term  in  the  integrand  suggests  that  the  main  contribution  to 
the  dose  at  point  D*  comes  fromiegions  where 


In  this  region,  all  fvinctions  except  are  slowly  varying  functions  of 

X.  We  therefore  can  approximate  them  by  their  value  at  x  =•  0,  ^(x)  = 


will  be  approximated  by  f(0)  + 


Hence 


X 

I 


Using  this  simplification  the  integration  over  x  can  be  easily  carried  out. 
We  obtain: 
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!  t 

l; 


The  integration  is  carried  out  over  the  area  which  can  be  seen  from  the  de¬ 
tector  at  D'.  ¥e  therefore  get  for  the  total  dose  D  at  T3' 


For  a  detector  placed  at  a  point  which  cannot  see  the  source,  the  direct  term 
is  eliminated.  and  the  relevant  areas  are  different  in  this  case,  however 
in  every  other  respect  the  approach  is  identical.  The  integral  over  the  sur¬ 
face  area  was  carried  out  numerically.  The  methods  of  the  numerical  integration 
will  be  explained  in  a  later  section. 


1 
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2,  Build-up  Factor  Correction  to  Single  Scattering 

We  shall  now  see  how  the  bulld-up  factor  would  effect  our 
single  scattering  results.  Ihe  bulld-up  factor  Is  defined  as  the  ratio 
between  the  amount  of  radiation  actually  arriving  at  a  certain  point  to  that 
wiilch  would  be  expected  should  there  be  no  scattering.  Scattering  builds  up 
radiation  at  a  point  because  scattering  is  not  a  removal  process.  The  amount 
of  radiation  arriving  at  the  scattering  volume  dV  Is 
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where  is  the  build-up  factor  and  is  effective  only  in  the  scattering 

material,  no-o  in  air.  Ihe  dose  at  the  detector  will  be  increased  again  while 
going  through  the  scattering  material  and  the  scattered  radiation  at  D*  will  be 


D  -  ^^2!* 


— p2t,c/i 


(21) 


The  integration  is  over  all  scattering  points  which  can  see  the  detector.  The 
functions  and  B2  are  not  identical  since  they  include  constai ts  which  depend 
on  the  energy  of  the  scattered  radiation  as  well  as  the  angle  of  incidence. 

In  the  case  of  gamma  rays  the  build-up  factor  can  well  be  approximated  by 


B(r(j)=i+/A(z:dr' 


(22) 
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Here  d  le  the  depth  the  radiation  penetrates  Into  the  scattering  material, 

A  Is  a  constant  and  n  an  Integer.  Both  A  and  n  depend  on  the  energy  and 
angle  of  Incidence  of  the  gamma  ray.  Introducing  (22)  into  (21)  we  obtain: 

(23) 


JV 


Dc  =  1-4 


TW^  K 


Making  the  same  approximation  as  was  done  for  single  scattering  and  elavuating 
the  integral  of  x  from  o  to  OO  we  obtain  for  the  scattered  doss 


The  contribution  from  the  direct  beam  is  added  to  (2[|).  The  angle  of  inci¬ 
dence  0  is  measured  from  the  vertical. 


-21- 


3.  Albedo  Approach 

We  shall  now  discuss  the  albedo  method  which  measures  the 
reflectivity  of  a  medium  with  respect  to  an  incident  flux  of  radiation.  It  will 
generally  depend  on  the  kind  of  radiation  under  consideration,  the  manner  in 
which  it  is  applied  (number,  energy  or  dose)  as  well  as  on  the  initial  condition 
(energy  and  angle  of  incidence).  We  shall  be  interested  in  dose  albedos  in  the 
case  of  gamma  rays,  and  in  number  albedos  in  the  case  of  neutrons.  In  the 
present  section  a  general  formula  is  developed  for  the  albedo  approach  and  for 
the  case  where  it  depends  on  both  the  initial  energy  as  well  as  on  the  angle 
of  incidence.  There  will  be  no  difficulty  in  applying  this  formula  to  neutrons 
for  a  constant  energy  albedo. 

Let  be  the  amount  of  a  certain  observable  quantity  (x 
will  indicate  the  specific  quantity  under  consideration)  falling  on  a  unit 
area  perpendicular  to  the  direction  of  the  beam.  Let  J*  (J2)  be  the  same 
quantity  emerging  per  unit  of  solid  angle.  Let  0  be  the  angle  of  incidence. 

Then  we  obtain: 


(25) 

Equation  (25)gives  the  distoibution  of  the  outgoing  radiation.  If  we  assume 
that  the  outgoing  distribution  is  isotropic  we  obtain 


f  {-^)  -  t^rr 


(26) 


Fig.  4.  Albedo  Geometry 


Figure  U  is  again  a  top  view  of  the  entranceway.  The 
source  is  located  at  5  while  the  detector  is  at  D‘.  Using  our  notation  we 
obtain 

r- 

WRf-  (: 

Here  S  measures  the  strength  of  the  source.  Using  (25),  (26),  and  (27),  w 
obtain  for  the  scattered  dose  at  D* 


The. manner  in  irtiich  this  integration  is  carried  out  will  be  indicated  in 
more  detail  in  a  later  section  treating  the  methods  of  integration.  There 
it  will  also  be  shown  how  to  take  account  of  multiple  reflection  whenever  it 
is  relevant. 

li.  The  Corner  Effect 

When  comparison  is  made  between  experiment  and  rigorous 
albedo  calculation  of  the  dose  distribution  in  a  duct,  Ihe  agreement  is  rather 
poor.  Albedo  calculations  result  in  lower  values  than  measured.  The  discre¬ 
pancy  is  worse  at  points  just  after  passing  the  bend  of  the  entranceway.  This 
suggests  that  there  exists  another  contribution,  which  so  far  has  been' 
neglected,  and  should  be  incorporated  in  any  calculation.  Considering  the 
experimental  facts  and  the  geometry  it  seems  quite  possible  that  direct 
transmission  through  the  corner  as  well  as  scattering  through  the  corner  is 
responsible  for  this  discrepancy.  We  therefore  are  interested  in  estimating 
the  possible  importance  of  the  corner  effect. 

Ihe  directly  transmitted  radiation  through  the  corner,  which 
decreases  exponentially,  is  small  even  for  the  first  position  after  the  bend. 

We  now  estimate  the  value  of  the  scattered  component. 

Unfortunately,  the  useful  albedo  concept  is  not  applicable 
directly  to  corner  scattering.  On  the  other  hand,  we  know  that  single  scat¬ 
tering  gives  a  considerable  underestimate.  We  now  introduce  a  fictitious  albedo 
to  determine  the  scattered  contribution  from  the  corner.  We  postulate  that  the 
ratio  between  contributions  from  two  different  sources  (wall  and  corner)  does 


-24- 


The  integration  is  over  the  surface  of  the  wall,  We  shall  now  perform  a 
similar  calculation  for  the  corner. 


Prom  Fig.  5.  we  can  see  that  this  contribution  Trill  be 


Here 

d,  —  P  ^ 


and  z  is  a  coordinate  perpendicular  to  the  drawing.  The  main  contribution 
v/ill  come  from  volume  elements  which  satisfy 
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Fig.  5.  CORNER  SINGLE  SCATTERING  GECMETRY 


This  allows  us  to  follow  the  approximation  developed  for 
single  scattering  from  a  wall  with  the  advantage  that  we  can  integrate  over 
two  coordinates,  ’.Te  obtain 

Cs 

Here 


We  are  left  with  the  integration  over  one  coordinate  only.  The  last  integral 
can  be  estimated  in  a  similar  way  to  that  of  equation  (30).  We  obtain 

Cs 

Using  (33)  and  (30)  we  obtain  for  a  solution  to  (29) 


^  (9)cornt^S^  I 

\  h^/i\ 


Z{0)  WON 


(3U) 
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Here  (0)  corner  is  the  average  aea tiering  angle  for  tha  corner  and  (6)  wall 
is  the  corresponding  quantity  for  the  wall.  How  this  correction  term  was 
uilllaed  will  be  shown  In  the  section  discussing  numerical  methods. 

B.  Calculated  Gamma  Hay  Albedo  Data 

The  fundamental  assumption  of  earlier  gamma  ray  streaming  theories 
was  that  the  albedo  Is  constant  with  respect  to  Incident  angle.  This 
assumption  has  been  proven  false  and  In  fact  the  angular  dependence  of  the 
albedo  Is  very  large.  In  obtaining  the  dose,  the  need  for  an  angular  dependent 
and  energy  dependent  albedo  la  deemed  entirely  necessary.  The  fundamental 
source  of  this  data  is  a  report  put  out  by  the  National  Bureau  of  Standards. 
Although  this  report  did  not  give  a  complete  table  for  variations  in  angle  and 
energy  it  did  cover  the  incident  angles  of  0*^,  30°,  60°,  and  90°  (grazing),  and 
energy  from  0. 2  to  2  mev.  The  variations  in  angle  and  energy  permit  fairly 
accurate  interpolation  and  extrapolation  for  any  angle  between  0  and  90°  and 
any  energy  between  0.1  and  3  mev.  With  regard  to  gamma  ray  exit  angle  varia¬ 
tion,  there  is  very  little  data  available.  Ihe  N.  B,  S.  Report  does  mention 
the  variation  of  the  dose  with  exit  angle,  but- the  data  is  sketchy  and  only 
general  conclusions  can  be  drawn.  The  other  reference  used  to  find  the 
suitable  albedo  is  by  J.  F,  Perkins Hiis  article  has  some  data  on  high 


Berger,  M.  and  D.  Raso.  "Backscattering  of  Gamma  Rays."  N.  B.  S.  Reoort 
No.  5982  (July  2h,  1958). 

Perkins,  J.  7.  "Monte  Carlo  Calculation  of  Gamma-Ray  Albedo  of  Concrete 
and  Aliiminum. "  Jour,  of  App.  Physics.  Vol.  l6.  No.  U,  (June  1955). 


tncrgy  albedos  (0,6  mev  to  6  mar).  Both  sets  of  data  were  developed  from  a 
Monte  Carlo  calculation  performed  on  a  computer.  N.  B.  S.  used  an  IBM  70h 
and  Perkins  used  an  19^  701.  The  N.  B.  S.  used  some  ^,000  case  histories, 
vdiereas  Perkins  used  only  2$0  case  histories  for  most  of  his  points.  Both 
methods  took  only  compton  scattering  and  absorption  into  account  in ’the  compu¬ 
tation  of  the  albedo.  Since  the  angular  dependence  of  the  emitted  gamma  ray 
requires  more  than  $,000  case  histories  in  order  to  obtain  sufficient  accuracy, 
the'  emitted  gamma  distribution  was  assumed  to  be  isotropic  in  our  calculations. 

In  particular,  the  albedo  depends  on  the  -incident  angle,  the  exit  angle,  and 
the  angle  between  the  incident  plane  and  the  exit  plane  as  illustrated  in  Fig.  6. 
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Figure  7  shows  plots  of  the  various  dose  albedo  data  used 


in  the  entranceway  problem. 

C.  Numerical  Techniques 

e 

1.  Duct  Equations  in  Rectangular  Geometry 
•  As  all  acperimental  work  is  done  in  rectangular  entrancerays  we 

present  the  equations  for  rectangular  geometry. 

The  expression  for  single  scattering  in  rectangular  geometry  is 


(x+s)%a +2^  >  2  -  (X.+  s)b^r 


.  -v  X  _p'  J - X - 3C  d' 


<lV=J 


So  the  scattered  dose  D  at  point  D’  (Fig.  8)  is 


r\  -AS  IM  P(6) .  « ? 

UTT^ 


0 
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Carrying  the  numerical  integration  for  then  integrating  over  x  at  x  -  0 

we  get 


irtiere  the  0,9  factor  appears  from  the  integration  which  will  allow  us  to  use 

2^22  22  2 
R  "  S  y  instead  of  ^  sec  'y  *  7  •  Doing  the  same  thing  with  the 

cylindrical  geometry  where  there  is  no  need  for  y  y  we  obtain 

(■f 

0 

■yilhen  the  build-up  factor  expression  B(  2.,  d)  *  1  +  A(21.d)^ 

"t  Xi 

is  included  in  the  single  scattering  formula,  the  final  expression  in  the 
rectangular  geometry  will  be: 


(37) 
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The  same  thing  is  done  for  the  ^vall  reflection  expression  (albedo  approach) 
where  we  get 


(38) 


. 

Here  oC=  -  to  account  for  the  multiple  reflection  from  the  four  walls. 

/"*  o4 

2.  Methods  of  Numerical  Integration 

The  expression  derived  above  can  be  numerically  integrated 
along  all  the  wall  surface  points  in  the  first  leg  to  find  the  reflected  contri¬ 
bution  to  any  point  on  the  centerline  of  the  leg  .  In  the  concrete  tunnel 
Ay  was  taken  to  be  one  foot.  and  are  calculated  at  every  point.  Also 
in  the  case  of  the  albedo  approach  and  2g(®)  and  P(0)  in  the  case 
of  the  single  scattering  are  de- ■  "mined  for  each  point.  In  ^  the  four  v;all 
scattering  surfaces  can  be  considered  symmetric  with  respect  to  a  dose  point 
at  the  centerline  in  spite  of  the  r  angle  dependence  in  our  geometry  which 
is  taken  care  of  by  the  factor  0,9.  This  expression  can  also  be  used  to  cal¬ 
culate  the  wall  reflection  contribution  to  the  centerline  dose  points  of  the 
second  leg.  Dose  points  in  the  leg  A  do  not  see  the  source  because  of  the 
right  angle  bend,  but  they  receive  a  reflected  contribution  from  all  wall 
scattering  sxirfaces  that  can  be  seen  by  the  dose  point  and  the  source  at  the 
same  time.  Rj^  and  Rg  are  calculated  from  each  of  these  scattering  surfaces 


to  the  source  and  tha  centerline  dose  point  specified.  It  should  be  borne 


in  mind  that  only  the  ceiling  and  the  floor  scattering  surfaces  at  the  comer 
are  symmetric  with  respect  to  the  centerline  dose  point. 

P(Q)  are  calculated  for  each  scattering  surface  also. 

In  addition  to  this  scattered  contribution,  dose  points 
in  the  second  leg  and  especially  those  near  the  corner  will  receive  a  con¬ 
tribution  from  the  radiation  that  falls  on  the  comer  and  scatters  in  their 
direction.  This  corner  effect  is  calculated  by  using  single  scattering  plus 
a  build-up  factor  as  described  previously. 

The  expression  for  the  corner  effect  correction  that 
should  be  added  to  the  albedo  contribution  is 


It  is  very  important  to  notice  that  the  reflected  contribu¬ 
tion  considered  here  from  the  points  that  can  be  seen  by  the  detector  and  the 
source  is  a  first  order  reflection  which  .contributes  directly  to  the  detector. 

A  second  order  reflected  contribution  can  be  considered  which  is  due  to  the  first 
order  reflected  radiation  being  scattered  by  the  second  lag  walls  to  the  detector. 
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fhia  la  done  by  considering  the  first  order  reflection  points  as  a  second 
order  source.  This  contribution  Is  not  going  to  be  important  In  the  ease  of 
gamma  rays  because  ve  are  measuring  dose  which  depends  on  energy  aid  the 
gamma  ray  energy  after  all  these  reflections  has  degraded  by  a  large  amount. 
This  contribution  will  be  of  extreme  importance  in  the  case  of  the  neutrons 
especially  when  we  know  that. the  detector  we  used  is  measuring  flux  per  unit 
time  and  the  albedo  we  used  for  this  reason  is  a  number  albedo.  The  first 
order  reflection  contribution  to  a  centerline  point  in  the  second  leg  is 


is 

37r^ 


and  the  flux  falling  on  a  unit  area  of  the  second  leg  wall  at  a  distance 
(Fig.  9.)  from  the  second  order  source  is 


and  the  contribution  from  this  point  to  the  detector  at  a  distance  R^  from  it 


kSjA  u'  _ 

/-uf  ZITIZ^ 


With  X  as  the  corner  effect  correction  factor  the  total  contribution  will  be 
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Cqo^JA 


The  express J.on  R2 


Cn<^dA 

■  u  "T^frtToir*  related  easily  to  the 

Z"”®' 


J 

integration  we  did  in  the  first  leg.  Here  the  number  albedo  was  chosen  to 
be  constant  for  different  experiments  done  in  concrete  and  found  to  be  0,675. 

3.  Sample  Calculations 

As  an  example  of  the  way  we  did  our  calculation  and  numeri¬ 
cal  integration  the  calculation  for  the  dose  at  position  No,  6  ig  demonstrated 
here.  As  it  is  shown  in  Fig.  9.,  this  point  is  12  feet  from  the  source.  The 
source  is  located  two  feet  from  the  duct  entrance.  The  source  used  is  a  3.67 
curie  Co^  gamma  ray  source  with  1,59  roentgen  per  hour  per  curie  at  one  yard. 
The  direct  dose  at  this  point  is 

1.59  X  (3/12)^  X  3.67  X  ^  »  6.032  mr/min. 

The  reflection  contribution  will  be  due  to  twice  the  area 
of  points  (1)  to  (3),  four  times  the  area  of  points  (Ii)  to  (lit),  three  times 
the  area  from  points  (l5)  to  (20)  and  the  area  frtsn  points  (21)  to  (28).  At 
point  8  for  example 

rJ  -  15.25  ft?  Rg  -  99.25  ft?  Cos  0  =  .768,o(=  .050, 

and  — =  .000026[i6  ft"^ 
l-0<  ^1  ^2 
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% 


0,002081i37 


but 


.1^^-  — —  7li2, 697  inr/min  •  ft^ 

^TT  Tf 


.  no-fj 

^TT  TT 


ieo^  71^2.697  X  .00208lili 


—  1. $h8  mr/min 


The  total  dose  D  at  position  No,  6  Is 


D  ”  6.032  +  l,51i8  ■  7.38  mr/min  . 

The  dose  without  the  roof  can  be  calculated  by  taking  the  roof  contribution 
which  is  due  to  points  (l3)  to  (20)  from  the  total  reflection  contribution 
which  makes  it 

71^2.697  X  0,0018U0  ■  1.3666  mr/min 
the  total  dose  =  6.0323  +  1. 338U  *  7,k99  mr/min 
The  same .thing  will  be  done  to  any  gamma  ray  source  with  only  one  variable 
whicli  is  the  albedo  value.  In  the  neutron  case  the  number  albedo  is  taken  as 
a  constant  as  we  have  no  energy  or  angle  dependence. 


IV.  BCPERniENTAL  PROGRAM 


The  experimental  part  of  this  program  places  major  emphasis  on  the 
measurement  of  the  attenuation  In  a  concrete  entranceway.  The  exact  size 
of  the  entranceway  Is  chosen  to  be  variable  so  as  to  apply  to  vent  ducts  as 
well  as  personnel  entraneeways.  With  these  ideas  in  mind,  the  entranceway  is 
constmcted  out  of  double- interlocking  concrete  blocks.  While  these  blocks 
were  being  obtained,  a  small  duct  with  a  right  angle  bend  was  constructed  out 
of  lead  to  check  theory  and  experiment.  More  work  on  lead  is  planned, 

A,  Description  of  Experimental  Entraneeways 

1.  Lead  Duct 

The  lead  duct  was  constructed  of  8”'  x  U*’  x  2”  lead  bricks.  The 
duct  has  an  8"  x  8”  cross  section  and  each  leg  is  16  ”  long,  A  photograph  of 

the  duct  is  shoivn  in  Fig,  10.  where  the  roof  has  been  removed  in  order  to  see 

the  actual  construction  in  more  detail.  The  walls,  floor,  and  roof  have  a 
minimum  thickness  of  four  inches  of  lead.  Note  ttiat  the  right  angle  corner 
has  been  filled  with  more  lead  brick  in  order  to  insure  no  direct  leakage 
through  walls. 

2,  Concrete  Entranceway 

The  concrete  blocks  are  designed  to  interlock  in  two  directions 
so  as  to  prevent  leakage  through  cracks.  The  blocks  ivere  of  various  sizes  and 
shapes,  and  range  from  an  8  foot  by  U  foot  by  1  foot  block  to  a  1  foot  cube. 

The  minimum  entranceway  size  possible  with  these  blocks  is  a  1  foot  by  1  foot 


cross  section  and  1  foot  length  to  a  6  foot  cross  section  and  2  foot  length. 
The  maximum  leg  length  is  about  30  feet  with  intermediate  cross  sections 
of  two  and  four  feet.  The  cross  section  does  not  have  to  be  square,  however 
•  this  simplifies  the  calculations.  The  drawings  in  Fig.  11,  show  the  general 
construction  of  the  concrete  blocks.  The  blocks  that  are  required  are: 

1)  standard  block,  2)  starting  block,  3)  ending  block,  U)  top  block,  5)  bottom 
block,  6)  start-top  block,  7)  start-bottom  block,  8)  end-top  block,  and 
9)  end-bottom  block.  All  blocks  are  derived  from  the  standard  block  by 
passing  a  plane  through  the  standard -block  perpendicular  to  the  flat  outside 
surface  as  is  noted  in  Fig,  11.  All  parts  of  the  wall  interlock  to  prevent 
leakage;  however,  the  top  row  and  the  bottom  row  of  blocks  are  flat  to  accom¬ 
modate  the  roof  and  floor.  The  parts  of  the  r oof  can  be  seen  to  interlock  in 
Fig.  12.  The  roof  itself  interlocks  but  not  with  the  wall.  Figure  13  shows 
the  concrete  entranceway  partially  assembled  with  the  reactor  thermal  column 
door  part  way  out.  The  entranceway  was  constructed  so  that  the  thermal  column 
door  can  pass  down  the  entranceway  and  the  reactor  used  as  a  neutron  (or 
gamma)  source.  Figure  lU  shows  the  almost  completed  entranceway  with  the 
right  angle  section  of  the'  roof  removed.  The  right  angle  roof  section  that  is 
shown  removed  is  the  physical  setup  used  when  the  data  in  this  report  refers 
to  "with  roof  removed, "  Figure  l5  shows  the  complete  entranceway  as  it  is 
assembled  for  the  measurements.  The  composition  of  the  concrete  is  as  follows 
To  one  bag  of  portland  cement,  205  pounds  of  sand  torpedo  l/dry,  2i;0  pounds  of 
3A  inch  stone  chips,  and  5  to  5.5  gal.  of  water.  i;,000  psi  compressive 
strength  resulted. 


B,  Gamma  Ray  Measurements 

1.  De tec tors 

Gamma  ray  dose  measurements  are  made  using  the  tandsverk 
Roentgen  Ueter  Model  This  set  of  intercalibrated  ionization  chambers 

is  capable  of  measuring  gamma  dose  from  a  Abaction  of  a  milliroentgen  to 
one  thousand  roentgens.  The  chambers  have  been  intercalibrated  by  the 
manufacturer  according  to  techniques  laid  dowi  by  the  Bureau  of  Standards, 

The  manufacturer  claims  a  precision  of  -  2  °/o,  and  provides  calibration 
data  for  correction  purposes. 

The  chambers  are  mounted  on  ring  stands  and  all  measurements 
are  made  on  the  geometric  centerline  of  the  entranceway.  Sufficient  care 
is  taken  in  positioning  the  detectors  and  correcting  for  calibration  and  back¬ 
ground  so  that  the  reported  dose  rates  are  accurate  to  better  than  -  10  °/o, 

2,  Sources 

•  The  gamma  sources  used  were  point  sources  of  Co^^  (3,67  curies) 
137 

and  Cs  (1,52  curies).  The  exact  calibration  procedure  is  described  later. 
The  Co^*^  sowce  was  held  in  a  cylindrical  brass  container  (1  in.  o.d.  x  8  in. 
long)  which  resided  in  a  lead  pig  when  the  source  was  not  being  used.  The 
cylinder  cover  was  attached  to  heavy  string  which  then  passed  through  a 
series  of  eyelets  fastened  to  the  roof  of  the  entranceway.  The  source  was 
moved  to  the  geometric  center  of  the  entranceway  by  operating  the  string  at  a 
safe  distance  from  the  source.  The  source  was  accurately  positioned  by 
sighting  on  a  second  string  mounted  horizontally  across  the  entranceway  at  a 
height  midway  between  the  floor  and  ceiling. 
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137 

The  Cs  source  Is  encapsulated  on  the  end  of  a  brass  rod 
of  3/16  in.  diameter  and  ll;  in.  length.  This  rod  is  fitted  with  heavy  string 

60 

and  is  positioned  in  the  same  way  as  the  Co  source, 

3.  Results 

Gamma  measurements  of  dose  rate  in  lead  and  concrete  entranceways 
are  summarized  in  the  tables  and  figures  to  follow.  The  data  is  presented  in 
such  a  way  as  to  be  more  or  less  self  explanatory.  Position  numbers  are  used 
as  a  handy  means  of  identifying  #iere  a  measurement  was  made.  (Figure  I6. ) 
Distance  along  the  centerline  measures  the  straight  line  distance  between 
source  and  detector.  ■  Use  of  the  dose  attenuation  factor  always  defines  the 
ratio  of  dose  rate  at  position  No.  6  (center  of  right  angle  bend)  to  the  dose 
rate  at  position  No,  1  (center  of  the  exit  plant).  Note  that  this  attenuation 
ratio  is  meaningless  without  the  ratio  .  being  specified. 

The  data  of  Table  II  is  for  measurements  with  a  point  source 

137 

of  Cs  gamma  rays  in  an  8  x  8  square  inch  cross  sectional  lead  duct  with 
U  inch  thick  walls  and  one  right  angle  bend  (see  Fig,  10), 

Table  III  summarizes  the  measurements  with  a  point  source  of 
Co^^  gamma  rays  in  the  same  lead  duct. 
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Fig.  16.  GAMMA  RAY  MEASUREMENT  POSITIONS 


-51- 


TABLE  tl 


POINT  SOURCE*  -  CAMMA  RAY  MEASUREMENTS 

IN  8'*  It  8"  LEAD  DUCT  WITH  4"  WALLS 


Position 

Number 

Centerline 
di stance 
from  source 
(inches) 

A 

(inches)  ' 

. 

(inches) 

Gamma 

Dose 

Rate 

mr/min 

Dose 

attenuation 

factor 

6 

20 

16 

0 

36.  58 

•  •  •  •• 

1 

40 

16 

16 

0.0537 

680 

L52  curies  source. 


TABLE  ni 

POINT  SOURCE*^-  GAMMA  RAY  MEASUREMENTS 

IN  8"x8"  LEAD  DUCT  WITH  4"  WALLS 


Position 

Number 

Centerline 
distance 
from  source 
(inches) 

(inches) 

(inches) 

Gamma 
Dose  ' 
Rate 
mr/min 

Dose 

attenuation 
fa  ctor 

6 

20 

16 

0 

331.  3 

•  • 

1 

40 

16 

16 

0.  766 

432 

3.  67  curie  source 


!,[ 

i 

L 

1. 

L 

.1. 

L 

[ 

L 

L 


The  dose  attenuation  factor  is  observed  to  be  some  36  ^/o 

smaller  for  Co^.  This  appears  explainable  on  the  basis  that  at  Co^^energies 

137 

the  Compton  scattering  is  the  dominating  process  vhile  for  Cs  energy  the 

photoelectric- absorption  process  Shares  about  equally ‘vdth  compton  scattering. 

Table  IV  summarizes  the  measured  data  for  a  point  source  of 
137 

Cs  gamma' rays  in  a  6  by  6  foot  square  cross  sectional  concrete  entrance- 
way  with  one  (1)  foot  thick  walls  (see  Figs,  10  through  lU). 

Table  V  summarizes  the  measured  data  for  a  point  source  of 
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Co  gaioma  rays  in  a  6  by  6  foot  square  cross  sectional  concrete  entrance¬ 
way  with  one  (1)  foot  thick  walls. 

Referring  to  Tables  IV  and  V  we  note  in  comparing  the  dose 
attenuation  ratios  in  concrete  we  have 


Source 

-4 

(feet) 

^2 

(feet) 

Bose  attenuation 
ratio 

Cs"37 

9 

16 

266 

Co6°  • 

9 

16 

h2h 

where  from  Tables  II  and  III  in  lead  we  found 


Source 

A 

( inches ) 

A 

(inches) 

Dose  attenuation 
ratio 

16 

16 

680  •  ■ 

o 

o 

c^ 

o 

16 

16 

1432 
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TABLE  IV 
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POINT  SOURCE  ,  Cs  GAMMA  RAY  DOSE  MEASUREMENTS 
IN  6'  X  6'  CONCRETE  ENTRANCEWAY  WITH  1'  WALLS 


L 

L 

Position 

Number 

Centerline 
distance 
from  source 
(feet) 

(feet) 

4 

(feet) 

Dose 

Rate 

mr/min 

Dosb 

attenxiation 

factor 

L 

8 

5 

9 

0 

4.  06 

____ 

L 

6 

12 

9 

0 

0.  818 

1 

4 

19 

9 

4 

0. 0401 

20.4 

L 

3 

23 

9 

8 

0. 01193 

68.5 

1 

1 

31 

9 

16 

0. 00307 

266 

* 

1.  52  curie  source 
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TABLE  V 
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POINT  SOURCE  .  Co  GAMMA  RAY  DOSE  MEASUREMENTS 


IN  6'  X  6*  CONCRETE  ENTRANCEWAY  WITH  V  WALLS 


Position 

Number 

Centerline 
distance 
from  source 
(feet) 

(feet) 

(feet) 

Dose 

Rate 

mr/min 

.  Dose 

attenuation 

factor 

7 

5 

5 

0  • 

33.96 

6 

8 

5 

0 

12.  19 

5 

11 

5 

0 

11.  16 

4 

15 

5 

4 

0.  6223 

19.5 

3 

19 

5 

8 

0. 1965 

62.  1 

2 

23 

5 

12 

0.  0793 

154 

•  1 

27 

5 

16 

0.  0403 

302 

7 

7 

7 

0 

17.  76 

6 

10 

7 

0 

7.  713 

5 

13 

7 

0 

7.  759 

4 

17 

7 

4 

0. 3660 

21.  1 

3 

21 

7 

8 

0.  1216 

63.4 

2 

25 

7 

12 

0. 04557 

169 

1 

29 

7 

16 

0.  02311 

334 

.  8 

■4 

9 

0 

53.4 

7 

9  . 

9 

0 

9.  21 

6 

12 

9 

0 

6.62 

,  5 

15 

9 

P 

5.  77 

4 

.  19 

9 

4 

0.  256 

25.  9 

3 

23 

9 

8 

0. 0761 

87.  0 

2 

27 

9 

12 

0.  0298 

222 

1 

31 

9 

16 

0.  0156 

424 

3.  67  curie  source. 
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We  thus  observe  the  interesting  fact  that  the  ratios  have 

reversed,  Co^  attenuation  being  larger  than  in  concrete  but  smaller 

in  lead.  If  we  observe  a*  plot  of  the  compton  and  photoelectric  cross  section 

for  concrete  we  find  that  both  for  and  Co^  the  compton  process  alone 
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acts, with  Co  having  a  small  scattering  cross  section  and  hence  a  greater 
attenuation  in  the  entranceway. 

The  lower  third  of  the  data  in  Table  V  is  plotted  as  Fig.  17. 

The  solid  line  represents  a  rigorous  albedo  calculation  including  corner  effect. 
The  dashed  line  shows  the  fraction  of  the  total  which  is  due  to  the  scattered 
■  component.  The  circled  points  are  measured.  In  general  the  agreement 
between  calculation  and  measiarement  is  rather  good.  Once  the  right  angle 
bend  is  passed  the  difference  between  scattered  and  total  is  due  to  the 
corner  effect  contribution.  It  is  not  clear  why  the  calculations  show  a  trend 
to  turn  upward  at  position  No.  1  (31  feet  from  source),  resulting  in  a  dis¬ 
agreement  with  measurement  of  approximately  U6  °/o. 

Table  VI  summarizes  the  calculated  data  for  the  6'  x  6'  concrete 
60 

entranceway  with  a  Co  point  source.  These  data  may  be  compared  with  the  data 
in  the  lower  third  of  Table  V  ( -9  feet).  Note,  that  the  calculated  (376) 
vs  the  measured  (i;2ii)  dose  attenuation  ratios  are  different  by  approximately 
11  °/o.  Reference  to  Fig.  17  shows  why.  At  the  31  foot  position  (No.  1)  the 
discrepancy  between  calculation  and  measurement  is  h6  °/o.  Note  that  no 
normalizing  of  data  is  done.  Dose  calculations  are  absolute. 
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3.  67  curies 
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The  calculated  gamma  dose  rate  has  three  components;  direct, 

Xor  those  positions  which  see  the  source,  scattered,  at  all  positions  and 
comer  effect  cintrlbutlons  which  decreases  as  leg  lengthens.  Note  that 
at  position  No.  U  (^2  •  U  feet)  the  corner  effect  contribution  (0.0596 

mr/min)  is  25  °/o  of  the  total  (0.2382  mr/min). 
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Referring  to  Fig.  18.,  for  Cs  ,  at  position  No.  1  almost  no 
discrepancy  exists  between  measurement  and  calculation.  The  measured  data  is 
taken  from  Table  IV  and  the  calculated  results  are  by  rigorous  albedo  theory 
including  corner  effect.  Note,  that  in  both  Fig.  17  and  Fig.  18,  the  calcu¬ 
lated  direct  dose  in  leg  is  higher  tlian  measured.  This  may  be  due  to  an 
error  in  source  strength  calibration.  Tiie  calibration  was  done  inside  the 
reactor  room  iflhere  scattering  off  walls  would  give  a  higher  indicated  source 
strength. 

In  order  to  verify  the  expected  effect  of  a  radiation  trap  the 
entire  roof  section  over  the  right  angle  bend  is  removed  as  shown  in  Fig.  IJ4. 
The  experimental  data  is  presented  in  Table  VII  aid  plotted  in  Fig.  19. 

The  effect  of  the  trap  may  be  seen  by  comparing  the  appropriate 
data  of  Table  V  and  VII,  From  Table  V,  without  a  trap  ,  the  dose  attenuation 
ratio  is  Ii2li  while  with  the  trap  the  ratio  is  5l7.  Thereof  section  con¬ 
tributes  about  equally  with  the  floor  and  side  wall  to  the  scattered  radiation 
because  the  angle  of  the  incident  radiation  is  the  same  and  near  90°.  For 
the  rear  scattering  surface  the  incidence  is  perpendicular  (zero  degrees)  and 
contributes  less  to  the  scattered  radiation  than  either  of  the  other  three 
surfaces. 
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Fig.  18  Comparison  of  Calculated 

and  Measured  C8-137  Gamma  Ray 
Dose  Distribution  in  6'  x  6' 
Concrete  Sntranceway 

point  source  -  1.52  curies 


Centerline  Distance  from  Source  (ft) 


.  TAB1.E  VII 
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POINT  SOURCE  -  Co  GAMMA  RAT  DOSE  MEASUREMENTS 
IN  6'  X  6'  CONCRETE  ENTRANCEWAY  WITH  1‘  WALLS  AND  WITH  TRAP 


Position 

Number 

Centerline 
distance 
from  source 
(feet) 

1 

(feet) 

2 

(feet) 

Dose 
Rate 
mr /min 

Dose 

attenuation 

factor 

8 

4 

9 

0 

51.9 

7 

9 

9 

0 

10.  38 

6 

12 

9 

0 

6.61 

4 

19 

9 

4 

0.  212 

31.  1 

3 

23 

9 

8 

0.  0574 

115 

1 

31 

9 

16 

0.0128 

517 

67  curie  source 
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Ihis  lack  of  equal  contribution  to  the  scattered  radiation 
going  down  leg  ^2  is  more  evidence  of  the  angular  dependency  of  the 
albedo.  Comparison  of  Figs.  17  and  1?  indicates  the  overall  effect  of  the 
trap.  The  results  of  the  trap  data  necessitated  reconsideration  of  earlier 
ideas  in  #iich  it  was  planned  to  p-lace  various  thicknesses  of  lead  on  the 
wall  immediately  at  the  bend.  It  was  argued  that  a  small  thickness  of  lead 
sheet  on  the  four  walls  would  greatly  increase  the  attenuation  of  the  right 
angle  duct  because  of  the  much  lower  lead  albedo.  The  lead  would  also 
reduce  the  corner  effect  materially.  This  is  planned  for  future  measurements. 

C.  Neutron  Number  Albedo  Measurements 

In  order  to  analytically  describe  neutron  number  flux  and  neutron 
number  attenuation  ratios  in  the  6*  x  6'  concrete  entranceway,  it  is  desirable 
to  again  resort  to  albedo  theory  because  of  the  success  this  receipt  provides 
in  describing  theganma  dose  distribution.  Unfortunately,  only  the  neutron 
dose  albedos  of  Strickler  and  Hungerford  are  available,  and  these  are  in  rather 
serious  disagreement.  (See  section  II-B). 

To  obtain  a  neutron  number  albedo  ,th8  following  two  measurements 
are  made;  A  PuBe  source  and  a  long  BF^  counter  are  mounted  on  12  foot  poles 
separated  by  six  feet  out  in  a  field.  An  average  of  twenty  readings' reatriWsS’’--  - 
in  a  count  rate  of  counts  per  minute.  Air  scattering  is  ignored  for  the 
6  foot  separation  of  source  and  detector. 
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The  same  source  and  detector  are  now  positioned  6  feet  apart 
and  on  the  axis  of  the  leg  of  the  concrete  entrancew^r.  The  position  of 

I 

the  counter  is  vertical  in  both  measurements.  We  obtain  a  count  rate  of  li58 

counts  per  minute.  As  we  have  the  same  direct  flux  as  outside  the  difference 

h^6  -  •  liOli  is  due  to  neutrons  scattered  from  the  walls.  We  now  must 

relate  these  two  measurements  to  the  neutron  number  albedo. 

Let  0  be  the  direct  flux  and  0  the  scattered  flux.  The  ratio 
0  s 

between  the  in  duct  and  free  air  measurements  is 


r-i 


(39) 


If  the  source  and  detector  is  identical  in  both  measurements  and  if  the  counter 
is  insensitive  to  the  energy  degradation  of  the  scattered  flux,  (assumed)  is 
the  count  rate  ratio  =  S.IjB.  From  previous  discussion  and  albedo  theory 


(b  - 


(Uo) 


where  is  the  source-detector  separation  and  is  the  source  strength  in 


neutrons  per  second. 


A 


(Ui) 
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which  is  from  equation  28  except  that  the  albedo  is  here  assumed  to  be 
independent  of  energy  and  angle.  is  the  distance  from  source  to  each 
wall  incremental  scattering  area  and  Rg  is  distance  from  wall  scattering 
area  to  detector. 

Forming  the  ratio  of  (UO)  and  (Ul)  we  obtain 


S 


The  integration  is  performed  numerically  as  described  in  the  sample  calculation 
of  section  III-C^3  with  the  results 


Hence,  from  equation  (39) 


ot=  0-675' 


This  value  for  the  neutron  number  albedo  is  used  in  all  calculations  to  follow. 
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D.  Neutron  Number  Flux  Measurements 


1.  Detectors 

Ideally,  as  for  the  gamma  rays,  a  measurement  of  neutron 
dose  distribution  in  an  entranceway  is  desired.  To  our  knowledge  no  accurate 
tissue  equivalent  neutron  dose  counter  is  commercially  available.  As  this 
program  lacked  sufficient  funds  to  develop  and  build  such  a  counter  an  attempt 

V. 

was  made  to  borrow  suitable  insti-uments  from  ANL  and  ORNL  without  success. 

We  thus  resorted  to  the  long  BF^  counter,  which  consists  of  a  BF^  counter 
surrounded  by  one  inch  of  a  paraffin  with  an  outer  cover  of  cadmium.  Neutrons 
with  energies  below  the  cadmium  cutoff  of  l.iiU  ev  are  prevented  from  being 
counted. 

Neutrons  of  energy  greater  than  l.Iiii  ev  in  energy  pass 
through  the  cadmium,  are  moderated  in  the  paraffin  and  detected  by  the  n,e<  ' 
reaction  in  the  EF^  counter.  The  overall  detector  assembly  measures  eight 
inches  in  length  by  three  inches  in  diameter.  The  active  volume  of  the  BF^ 
tube  is  U-l/2  inches  in  length  by  1  inch  in  diameter.  While  the  energy 
response  is  unknown, such  long  counters  characteristically  have  a  somewhat 
uniform  sensitivity  to  neutrons  from  a  few  raev  to  the  resonance  region. 
Suitable  amplifier  bias  removes  the  normally  small  gamma  sensitivity. 

It  should  be  mentioned  that  use  of  a  5  curie  PuBe  source 
results  in  fluxes  in  the  6'  x  6'  concrete  entranceway  which  are  too  small  for 


1 1 


.1'  .  i 

foil  threahold  detectora  or  tissue  equivalent  ionization  chambers.  Scin¬ 
tillation  crystals  fpr  accurate  neutron  dose  rasasurements  are  lonavailable 
commercially.  Certain  mixtures  of  ZnS  and  Incite  are  known  to  produce 
tissue  equivalent  doses  but  are  beyond  the  scope  of  this  program  to  date. 

The  proton  recoil  counter  developed  by  Hurst  and  associates 
at  ORNL  are  now  available  commercially  and  will  be  used  for  future  neutron 
dose  measurements. 

2.  Sources 

An  investigation  of  possible  neutron  sources  for  use  in 
entranceway  measurements  resulted  in  the  selection  of  Pu-Be  as  a  practical 
source.  Consideration  was  given  to  the  use  of  fission  plates  to  produce  both 
a  fission  gamma  and  neutron  spectra  but  were  discarded  because  of  the 
inordinate  amount  of  time  required  for  fabrication.  Design  of  a  suitable 
plate  for  future  use  has  been  completed.  This  plate  requires  cooling  and  will 
produce  approximately  10^*^  neutrons/sec. 

The  PuBe  spectrum  emits  neutrons  ranging  in  energy  fromrv^ 
0.1  to  12  mev  with  an  average  energy  of  approximately  U.5  mev.  This  is  above 
the  fission  average  of  approximately  2  mev  and  well  below  the  ll;  mev  fusion 
energy.  As  variation  in  dose  rate  with  neutron  energy  is  not  large  in  the  1  to 
10  mev  range, the  source  appears  to  be  a  suitable  substitute  for  both  fission 
and  fusion.. 


I 

I 


A  5  curi«  PuB#  source  vraat  procured  from  ANL  on  loan  and  is 
convertible  to  a  point  thermal  source  by  emersion  in  a  12  inch  paraffin 


cylinder, 

3,  Results 

The  small  lead  ducts  were  not  of  practical  size  for  neutron 
measurements.  Vary  small  sources  are  required  and  thick  walls  to  reduce 
to  zero,  fast  neutron  penetration  of  walls.  All  neutrons  measurements  are 
made  in  the  6*  by  6*  concrete  entranceway.  Because  of  its  very  low  gamma 
emission  and  average  neutron  energy  in  the  range  of  interest,  PuBe  is  chosen 
as  the  source.  Only  point  source  measurements  have  been  made  to  date  through 
a  number  of  other  sources  are  possible  in  future  work. 

PuBe  point  source  neutron  data  is  summarized  in  Table  VIII 
and  plotted  in  Figure  21.  For  convenience,  the  positions  in  which  measure¬ 
ments  are  made  may  be  located  on  Figure  20.  for  all  neutron  measurements. 

Referring  again  to  Fig.  21,  the  solid  line  is  for  the  calculated 
neutron  number  flux  distribution  in  the  6'  x  6*  concrete  entranceway.  The  cir¬ 
cled  points,  connected  by  a  dotted  line,  are  for  measured  points.  The 
calculated  distribution  is  determined  by  albedo  theory  including  multiple 
reflections  and  corner  effect. 

As  with  the  gamma  ray  measurements,  neutron  measurements 
are  made  at  the  indicated  positions,  with  the  entire  enti^apceway,  i. e.. 
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6.  Sn/cm  -sec  =  1  mrem/hr. 


Fie.  20.  NEUTRON  NIR'BER  FLUX  I.EASUREMENT  POSITIONS 
IN  6»  X  6'  CONCRETE  EN'IRANCEMY 
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^2  -  16"  feet  always  present.  Hence,  Table  VIII  indicatee  an  X’g  length  as 
If  the  measurement  is  made  with  that  positten  at  the  end  of  The  presence 

of  more  wall  scattering  area  behind  the  measurement  position  certainly  adds 
some  s^jiall  contribution.  The  conversion  factors  indicated  at  the  bottom 
of  Table  7III  are  derived  in  a  later  section  of  this  report. 

Table  IX  presents  the  calculated  data  plotted  in  Fig.  21. 
Comparison  of  calculated  vs  measured  neutron  number  attenuation  ratio  is 
seen  from  Tables  VIII  and  IX  (position  No,  1)  to  be  U3.7  calculated  and  26.1 
measured  for  an  of  16  feet  with  “  3.  We  suspect  the  measurement 

for  a  number  of  reasons.  Assuming  the  number  fliix  falls  of  as  we  have  measured 
and  calculated,  it  is  important  to  emphasize  that  the  dose  distribution  would 
fall  of  considerably  faster.  The  neutrons  at  the  exit  end  of  the  entranceway 
are  degraded  in  energy  and  therefore  deliver  less  dose.  The  small  number 
attenuation  ratios  are  never- the- less  worthy  of  some  concern  for  shelter 
designers  and  future  experiments  should  measure  a  neutron  dose  distribution, 
tissue  equivalent  if  possible. 

The  effect  of  a  trap  is  again  investigated  by  removing  the  roof 
section  over  the  right  angle  bend  in  the  6*  x  6’  concrete  entranceway.  The 
data  are  presented  in  Table  X  and  plotted  in  Fig,  22, 

A  comparison  of  the  neutron  number  attenuation  ratios  of  Tables 
VIII  and  X  indicates  almost  a  negligible  difference.  The  roof  scattering 
surface  and  hence  by  symmetry  the  floor  surface  makes  a  very  small  contribution 
to  the  number  flux  at  the  exit.  This  is  surprising  as  we  have  had  other  indi¬ 
cations  that  the  angular  dependence  of  the  neutron  number  albedo  is  small. 

This  evidence  is  to  the  contrary. 
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Source  =  7.  09  x  10  n/sec 


\ 

\ 


V 


Centerline  Distance  from  Source  (ft) 


TABLE  X 

MSZASURED  PuBe*  NEUTRON  NUMBER  FLUX  DISTRIBUTION 


IN  6'  X  6*  CONCRETE  ENTRANCEWAT  WITH  TRAP 


Position 

Centerline 

a 

C/min 

Flux** 

Number 

Number 

dietanee 

Attenuation 

from  source 

Ratio 

Cfeet} 

(feet) 

(feet) 

n/cm  >’sec 

I 

31 

9 

16 

la 

30 

9 

15 

.  54 

2.93 

20.  0 

2 

27 

9 

12 

77 

4.  17 

14.  0 

2a 

24 

9 

9 

120 

6.48 

9.  05 

3a 

21 

9 

6 

203 

11.0 

5.  33 

4a 

18 

9 

3 

356 

19.  3 

3.  04 

5 

15 

9 

0 

867 

47.  1 

1.  24 

6 

12 

9 

0 

1082 

58.  6 

7 

9 

9 

0 

1811 

98 

8a 

6 

9 

0  . 

3055 

165 

9a 

3 

9 

0 

5533 

300 

Source  =  7. 09  x  10  h/sec 
2 

16.  9n/cm  -sec  =  312C/min. 


75 


In  an  attempt  to  get  some  information  about  the  energy  spectra 
of  the  neutron  distribution  in  the  concrete  entranceyray ,  a  series  of  thermal 
neutron  measurements  are  made.  Table  XZ  presents  data  for  a  PuBa  point 
source  and  thermal  neutron  distribution  measurements.  The  data  is  plotted 
in  Fig.  23.  The  PuBe  source  is  immersed  in  paraffin  sufficiently  thick  to 
thermalize  all  neutrons.  With  this  thermal  source,  replacing  the  PuBe 
spectrum  source,  another  thermal  neutron  distribution  measurement  is  made. 
These  data  are  given  in  Table  XII  and  plotted  on  Fig.  23. 

From  Fig.  23;  two  points  of  interest  are  evident.  First,  the 
thermal  number  flux  attenuation  ratio  is  distrubingly  small  (  /\ylS  from 
Table  XII).  Recall  the  PuBe  number  flux  attenuation  ratio  was  also  small 
( 1^3. 7  from  Table  IX). 

Second,  the  thermal  flux  distribution  from  the  PuBe  source 
falls  off  with  a  very  similar  slope  to  the  thermal  flux  distribution  from  the 
theimial  source.  One  might  suspect  that  as  we  proceed  down  the  entrancsway 
the  number  of  thermals  would  increase  because  of  thermalization  effects. 

Little  more  can  be  learned  from  the  thermal  flux  data. 
Apparently  the  flux  thermalizes  significantly  early  in  leg  Hence,  ’ 

again  a  dose  distribution  is  needed  and  would  show,  based  on  these  results 
that  the  dose  attenuation  ratio  is  far  larger  than  the  measured  number  flux 
attenuation  ratio. 
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TABLE  XI 


i 

■I 

1 

,  MEASURED  THERMAL  NEUTRON  NUMBER  FLUX  DISTRIBUTION 


IN  6'  X  6*  CONCRETE  ENTRANCE  WAY  WITH  PuB*  SOURCE 


Position 

Number 

Centerline 
distance  to 

source 

(feet) 

/ 

1 

(feet) 

(feet) 

Counts/min 

Number 

attenuation 

ratio 

1 

31 

9 

16 

. 

■■  -  ■ 

la 

30 

9 

15 

145 

10.52 

2 

27 

9  • 

12 

235 

6.49 

2a 

24 

9 

9 

352 

4.33 

3a 

21 

9 

6 

485 

3.  15 

4a 

18 

9 

3 

820 

1.86 

5 

15 

9 

0 

1244 

^  m  m  m 

6 

12 

9 

0 

1527 

7 

9 

9 

0 

1910 

8a 

6 

9 

0 

2647 

9a 

3 

9 

0 

3390 

Source  =  7. 09  x  10  n/sec. 


-  77  - 


TABLE  Xn 


I 

1 
L 

i  MEASURED  THERMAL  NEUTRON  NUMBER  FLUX  DISTRIBUTION 

•  L 

IN  6»  X  6'  CONCRETE  ENTRANCEWAY  WITH  THERMAL  NEUTRON  SOURCE 


Position 

Number 

Centerline 
distance 
from  source 
(feet) 

A 

(feet) 

4 

(feet) 

Counts 

per 

minute 

Number 

attenuation 

ratio 

1 

31 

9 

16 

V  •  •  • 

la 

30 

9 

15 

85 

14.5 

2 

27 

9 

12 

147 

8.  38 

2a 

24 

9 

9 

232 

5.31 

3a 

21 

9 

6 

372 

3.  31 

4a 

18 

9 

3 

559 

2.  21 

5 

15 

9 

0 

990 

6 

12 

9 

0 

1233 

.... 

7 

9 

9 

0 

1826 

8a 

6 

9 

0 

2761 

9a 

3 

9 

0 

4497 

V.  SOURCE  CALIBRATION  PROCEIXJRE 

Aside  from  a  comparison  of  calculated  vs  measured  attenuation  ratios 
in  the  ducts  it  is  desirable  to  coitqsare  on  an  absolute  ba'is  the  calculated  vs 
measured  dose  distributions.  The  precision  of  the  Landsverk's  chambers  makes 
this  feasible  but  of  questionable  value  in  the  case  of  neutrons  measured  with 
a  long  BFj. 

A.  Gamma  Source 

Ihe  source  and  detectors  are  places  as  far  as  possible  from  any- 


scattering  surfaces.  Detector  readings  are  taken  at  various  distances  from 
the  source.  The  results  -are  plotted  in  Fig.  2h.  for  bo-bh  Cs^^"^  and  Co^*^ 


Fig.  2J4.  Gamma  Source  Strength  Calibration  Determination 


12 

Using  the  conversion  factor  1,32  H/hr  per  curie  at  1  meter 


60  /  137 

for  Co  and  0,356  R/hr  per  curie  at  1  meter  for  Cs  vre  find 


1 . 


Source 

Curies 

c.“ 

3,67 

0.W’ 

1.52 

dis/see. ) 

B,  Neutron  Sources 

The  PaBe  neutron  source  Is  calibrated  by  Mound  Laboratory  as 
to  its  neutron  emission  rate.  The  value  is 

7.09  X  10^  n/sec,  5  curies 

For  distances  in  which  air  scattering  is  negligible  ,the  conversion  from 
source  strength  to  mrem/hr  may  be  made  as  follows,  at  6  feett 


7.0?  X  10 


i*nR  Uw  (6  X  12  X  2,5U) 


-j—  ■  l6,?  n/cm  -sec 


From  the  N.  B.  S,  Handbook -No,  63,  the  flux  equivalent  to  1  mrem/hr,  F(E), 

at  all  energies  in  the  PuBe  spectrum  is  found  and  plotted.  An  average  value 
2 

of  P(n/cm  -sec  per  mrem/hr)  is  obtained  by  numerical  integration  of  the  PuBe 
energy  flux  distribution^^  and  the  PCE)  function. 


Radiological  Health  Handbook,  U.  S.  Sept,  of  Health,  Education  and 
Welfare  (Jan.  1957). 

Stewart,  L.  Phys.  Rev.  98  (1955). 
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,0(B)  P  (E)  dE 


0(E}  dE 


19  find 


2 

6.8  n/cm  -sec  ■  1  mrem/hr 


for  PtiBa. 


Hence,  at  a  distance  of  six  feet  in  air  m  have 

-  2,li9  mrem/hr. 


As  stated  previously,  in  the  neutron  albedo  measurement,  we  have  detector 
readings  at  six  feet.  The  results  are  summarised  below* 


Source 

(Curies) 


Yield 

(n/sec) 


Plux 

(6  feet) 
(n/cm^-sec) 


Dose  Detector  C/min 

(6  feet)  reading  per 

mr€Bi/hr  C/min.  rareo/hr 


7.09  X  10  16.90 


VI.  SDIMARY  AND  CONCL05I0NS 

The  material  reported  herein  describes  a  series  of  analytical  and 
experimental  efforts  directed  toward  the  determination  of  neutron  and  gamma 
ray  distributions  in  ducts  and  entranceways.  By  ground  rule  agreement, 
integral  transport  theory  and  Monte  Carlo  methods  are  excluded  from  use  on 
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this  program.  The  albedo  receipt  is  applied  but  with  considerable  modifi¬ 
cation  and  rigor.  The  agreement  of  tha  albedo  model  and  experiment  is 
rather  good  for  gamna  rays  and  rather  poor  for  neutrons. 

The  orientation  of  the  program  is  such  that  full  scale 
personnel  shelter  concrete  entranceways  are  emphasised.  A  6*  x  6'  concrete 
entranceway  with  1  foot  thick  walls  Is  constructed  with  a  single  right  angle 
bend  but  with  provision  for  a  second  bend.  Lead  duets  are  also  used. 

Gamma  dose  distributions  are  calculated  and  measured  in  both 
lead  and  concrete  ducts.  Agreanent  between  calculation  and  measurement  is 
good. 

Neutron  number  flux  distributions  are  calculated  and  measured 
with  considerably  less  agreement.  The  results  of  the  flux  measurements  clearly 
show  that  neutron  dose  measurements  must  be  made.  Such  measurements  are 
now  possible  and  with  considerable  accuracy  through  use  of  the  Stanrad  Neutron 
Dose  Detector. 

The  need  for  a  major  dose  albedo  experimental  effort  for  neutrons 
and  gamma  rays  is  of  paramount  Importance.  This  laboratory  Is  planning  such  a 
program  in  conjunction  with  further  entranceway  streaming  measurements.  Speci¬ 
fic  recommendations  for  additional  entranceway  studies  are  contained  in  a 
proposal  submitted  concurrently  with  this  report. 
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Bie  work  presented  in  this  appendix  is  a  sumraaiy  of  a  short  but 
rather  accurate  method  to  the  solution  of  the  square  cross  section  duct  problem. 
The  method  holds  well  over  a  range  of  length  to  width  ratios  of  1  or  2  up. 

Neither  the  exit  angle  or  the  angle  between  the  entrance  and  exit  planes  is  taken 
into  account  in  the  plots  of  albedos  on  page  32  because  of  the  lack  of  informa¬ 
tion,  This  approach  gives  the  angles  to  use  iirtien  this  information  becomes 
available.  The  approach  is  broken  down  into  three  contributions: 

1.  Albedo  Wall  Scatter 


A^g(area  of  wall  seen  by  souree)  =• 
A^gj(area  of  wall  seen  by  detector) 


(  -^2*  2 
- ^ 


irtiere: 

^2  refer  to  the  lengths  of  legs  1  and  2  respectively, 
H  is  the  height,  and 

^  2  ”  hsli’-'width  of  legs  1  and  2  respectively. 
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where  n  is  the  attenuation  coefficient  (neutron  or  gamma)  of  the  corner  material 


I 

[ 


and  0'  ia  the  angle  from  the  source  to  the  corner  relative  to  the  wall  of  the 
0 

duct}  and  is  actually  the  penetration  angle.  Since  the  dose  decreases  very 
rapidly  with  an  increase  in  x,  the  effective  added  width  of  the  tunnel  x 

os 

(Fig.  25.)  is  quite  small  and  0'  0.  Thus  the  total  dose,  iiuoar ted  to  the 

wall  (after  multiplying  by  the  similar  triangle  relation  (.■■ —  y  i  )  and 

A 

height  H  isi 


r 


Area  x  Dose 


*  2  J 


a  )  H  (Dose  at  point  y  e  ° 
''I  without  comer) 


-u  (x/sjn0^) 


r 


tl  X 

a. 


■) 


J 


/?1  2  J  sin  0 

( - - ^)  H  (D  )  - - 

yo 


dx 


Considering  D^,^  to  be  approximately  the  same  as  the  average  dose  to  the  total 

wall  surface  A  ,  we  have 
wss* 


A  .  (total  area  seen  by  source)  =  A  + 
stw  wss 


(sin  0(j)  H  (  '^1  2  S  g) 


A  ,  -  A 


stw 


wss 


P  sin  0  n 

I 


1 
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Similarly: 


area  seen  by  detector) 


(sin  H  (  ^2  *  2 


A  *  A 

dlsff  *irsd 


sin  0^ 


where  and  are  the  average  angles  to  the  corner  from  the  source  and 
detector  points  respectively,  and  and  [i  are  the  attenuation  coefficients  for 
the  incident  and  scattered  radiation  respectively. 

These  angles  0^  and  0^  are  closely  given  by: 


-1 

00  ■  tan  — )  and  0^  =  tan"^  (L,.^  )  ^ 


As  a  rule  of  thumb  -  the  scattered  radiation  has  an  energy  of  the  original 
scattered  radiation  after  going  through  a  90°  bend  (doiim  to  approximately  0.5 
mev  for  gamma  rays).  The  extra  corner  penetration  also  adds  areas  to  the  top 
and  bottom  thusly: 

Ag^^(total  area  of  the  bottom  seen  by  the  source)  =  (1^  ^2^ 

tan  0  ^  rs  » 

*  **  tr  '  Sj) 
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1 


tan  0  5 

tan  0^  +  ■  ■■— —  (1  + 


Tl’ 

«■ 


area  of  top  seen  by  detector)  ■  (U  Sx  Sj) 


5i‘4A) 

H  Xg 


\  Si  tan  0  Si 


The  average  reciprocal  squared  distance. to  these  jareas  (or  surfaces) 


is  given  by: 


Average  Reciprocal  Square^  Average  Reciprocal  Square^ 
Distance  from  Source  (l/Rg)  Distance  to  Detector  (l/R^) 


(  *  I  Si*  jHj  (  /j  *  *  Jh‘ 


r^TTsTT^sTF 


^stb.  ^dtt 


(  ♦  S2>"  *  t  “ 


t  -^2  *  Sl>''  *  C  H‘ 


The  angle  0  is  the  angle  from  the  source  (incident)  or  detector  (exit)  to  the 
average  point  of  the  surface  as  measured  from  a  perpendicular  to  that  surface. 


1 

I 
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The  final  form  for  the  albedo  scattered  dose  is: 


“d  ~  »d 

Atb  ^  Att>  '=°°  V, 


Where  k  is  the  flux  to  dose  conversion  factor.  The  form  vdth  only  incident 
angle  albedo,  such  as  those  derived  previously  would  bet 


“<o  At.  ®1  0<< V  .  At,“'»®2<<'*2> 


<*3tb  "  Att)  °°°  ®3 

2v  eJ 
s  d 


2.  Corner  Scattering  Toward  the  Detector 

Using  single  spattering  (the  differential  Klein-Nishina  cross 
section)  technique  (described  in  more  detail  in  this  report  under  "single 
scattering"  section),  the  amount  of  scattered  radiation  that  is  scattered 
tl'irough  an  average  angle  of  90®  -  -  0^  is: 


D 


ac 


N  Sk  sin  0  sin  0^ 
0  ox 

- 2 - 2 - 

hnr^  Vo^" 


K(90®  -  0^  -  0^,  E^)  P(90®  -  0^,  0^,  E^) 


where 
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rl-i'fi*  Si)^  +  3/2  Sl,±n 


and  n  in  cm  ,  H  in  cm,  and 


K  •  the  Klein-Nishina  differential  cross  section  for  scattering 
through  an  angle  of  90®  0^  -  0^.*  the  energy  of  the  incident  radiation  E^, 

in  cm“^  of  comer  material.  P  is  the  energy  degradation  factor  for  the  above 
angle  and  energy.  For  neutrons,  the  use  of  ^An  should  be  sufficient  for 
K,  and  P  almost  unity  for  heavy  elements. 

3,  Direct  Penetration 

For  small  ducts, direct  penetration  may  be  a  problem,  however  in 
general  it  is  not.  The  dose  would  be 


N  ke 
0 


3  - - ; - : 

Un  (  +  $■])  ^^^2  *  ^2^ 


r  is  equal  to  the  distance  traveled  through  the  material,  ^i^  =  total  absorption 
coefficient  ^or  dose  attenuation  coefficient  for  neutrons),  and  B^(^i^r)  is  the 


dose  build-up  factor. 


The  total  dose  is: 


D  -  Dg  +  Dg  +  Dp. 
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